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In the presence of a strong magnetic field and under condi- 
tions as realized in the crust and the superfluid core of neutron 
stars the Hall-drift dominates the field evolution. We show 
by a linear analysis that for a sufficiently strong large-scale 
background field depending at least quadratically on posi- 
tion in a plane conducting slab an instability occurs which 
rapidly generates small-scale fields. Their growth rates de- 
pend on the choice of the boundary conditions, increase with 
the background field strength and may reach 10 3 times the 
ohmic decay rate. The effect of that instability on the rota- 
tional and thermal evolution of neutron stars is discussed. 



In the presence of a magnetic field the electric conduc- 
tivity becomes a tensor and, what is more, two non-linear 
effects are introduced into Ohm's law: the Hall-drift and 
the ambipolar diffusion. However, if the conducting mat- 
ter consists of electrons and one sort of ions and no neu- 
tral particles take part in the transport processes the 
ambipolar diffusion is absent |jj . Such a situation is real- 
ized e.g. in crystallized crusts of neutron stars and/or in 
their cores if the neutrons are superfluid, but the protons 
are normal and the electrons may therefore collide with 
protons but effectively not with the neutrons. 

The effect of the Hall-drift on the magnetic field evo- 
lution of isolated neutron stars has been considered by 
a number of authors (see e.g. They discussed the 

redistribution of magnetic energy from an initially large- 
scale (e.g. dipolar) field into small-scale components due 
to the non-linear Hall-term. Though the Hall-drift is a 
non-dissipative process, the tendency to redistribute the 
magnetic energy into small scales may accelerate the field 
decay considerably. 

Indeed, when starting with a large-scale magnetic 
field the Hall-cascade derived in || will generate small- 
scale field components down to a scalelength l cr u , where 
the ohmic dissipation begins to dominate the Hall-drift. 
Considering numerically the evolution of a magnetic field 
in a sphere consisting of spherical harmonics up to a mul- 
tipolarity 1 = 5, Shalybkov & Urpin || concluded that 
the inclusion of higher harmonics will not influence the 
magnetic evolution. This conclusion is what we want to 
put in question. 

In some of the above-mentioned investigations numer- 
ical instabilities are reported if either too many harmon- 
ics were taken into account [JtJ or the initial field is too 
strong ||. Also, when considering the thermomagnetic 
field generation in the crust of young neutron stars || , 
where small-scale modes are the first ones to be excited 



numerical instabilities occurred caused exclusively by the 
Hall-drift. 

Here we want to show that all the observed instabili- 
ties are very likely in their essence not of numerical ori- 
gin but have physical reasons. In presuming that we felt 
strongly supported by the close analogy of the linearized 
field evolution equation including Hall-drift to the induc- 
tion equation including the so-called to x j— effect intro- 
duced by Radler |l0| . Within the framework of mean- 
field dynamo theory (see e.g. pJ|) he demonstrated the 
possible occurrence of magnetic instabilities in an electri- 
cally conducting fluid if a shear flow acts together with an 
electromotive force (e.m.f.) perpendicular to the current 
density j cx curlB. 

To prove the existence of a Hall-drift induced insta- 
bility we employ a simplified model: We assume spatial 
constancy of the conductive properties of the matter and 
show that under special conditions with respect to the 
(large-scale) background field strength and geometry an 
instability occurs which quickly transfers magnetic en- 
ergy from the background field to small-scale perturba- 
tions. We present the result of a linearized analysis which 
returns only growth rates and the spatial structure of the 
unstable field modes. Only a fully non-linear analysis is 
able to yield saturation values of the excited small-scale 
modes. 

The instability may work in different physical systems 
but is probably most efficient in modifying the field decay 
in compact astrophysical bodies. Then it will act only 
during an episode of the field decay, which unavoidably 
leads to a zero field. This episode, however, may have 
observable consequences. 

In the absence of motions and of ambipolar diffusion, 
the equations which govern the magnetic field are 



B = — c curl ( ( curl B + LL>BT e ( curl B x en )) 

\4ira K B eV 7(1) 

div B = , 

where c is the speed of light, a the electric conductivity 
caused by electrons, r e the electron relaxation time and 
wb = e\B\/m*c the electron Larmor frequency, with e 
being the elementary charge and m* the effective mass 
of an electron, eg is the unit vector in indirection. An 
estimate of the two terms on the r.h.s. of (Q) gives rise 
to the supposition that the Hall-drift becomes important 
only if LUBT e > 1. 

Using standard arguments one can immediately state 
that in the absence of currents at infinity the total energy 
of any solution of ([j]) is bound to decrease monotonically 
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to zero since the Hall-term oc curl-B x eg is unable to 
deliver energy (nor to consume it). 

For simplicity we assume the conductive properties of 
the matter to be constant in space and time, that is, we 
assume constant er and r e /m*. Thus, the induction equa- 
tion can be rewritten in dimensionless variables such that 
it no longer contains any parameter and the magnetic 
field evolution is solely determined by its initial configu- 
ration B(x, 0). For that purpose we normalize the spatial 
coordinates by a characteristic length L of the model (for 
a neutron star it could be, e.g., its radius or the thickness 
of its crust), the time by the ohmic decay time AiraL 2 /c 2 
and the magnetic field by — m* e c/eT e . The governing 
equations in these dimensionless variables read 



B = AB - curl( curl B x B) , divB = 



(2) 



where the differential operations have to be performed 
with respect to the now dimensionless spatial and time 
coordinates x,y,z and r, respectively. 

Stepping now into the search for instabilities we first 
have to define a proper reference state Bq. In order to 
avoid difficulties in defining the term "instability" and to 
facilitate the calculations we assume -Bo to be constant 
in time. Consequently we are forced to assume the ex- 
istence of an additional e.m.f. which prevents Bq from 
decaying. Although appearing to be very artificial, we 
find this measure to be legitimate as long as the results 
of the stability analysis are applied to real physical situ- 
ations obeying the constraint that the background field 
B is changing only slightly during the considered period 
of time. 

Linearization of (|^) about Bq yields 



Let us now specify the geometry of our model and the 
background field. We consider a slab which is infinitely 
extended both into the x- and y-directions but has a fi- 
nite thickness 2L in z-direction. The background field is 
assumed to be parallel to the surface of the slab point- 
ing, say, in x-direction and to depend on the z coordinate 
only, i.e. -Bo = f(z)e x . Guided by the conditions under 
which the above-mentioned magnetic instability [fiof may 
work, we conclude that f(z) has to be at least quadratic 
thereby ensuring that the first term in (|3[) is able to play 
the role of the shear flow, the second the role of the Co x j— 
term. Note, that by this choice curl -Bo x -Bo represents 
a gradient. Thus the unperturbed evolution of the back- 
ground field is not at all affected by the Hall-drift; in the 
absence of an e.m.f. it would decay purely ohmically! 

Further on we decompose a perturbation b into a 
poloidal and a toroidal component, b = b p + b t , which 
can be represented by scalar functions S and T, respec- 
tively, by virtue of the definitions 

b p = -cml(e z xWS) , b t = -e z xVT , (4) 

ensuring div b — for arbitrary S, T. 

For the sake of simplicity we will confine ourselves to 
the study of plane wave solutions with respect to the x- 
and y-directions, thus making the ansatz 



j (x, t) = | & \ (z) cxp (ikx + pr) , ( ) 



where k = (k x ,k y ), x = (x,y) and p is a complex time 
increment. It guarantees as well the uniqueness of the 
poloidal-toroidal decomposition since from A(S,T) = 



b = Ab - curl( curl B x b + curl 6 x B ) , div& = it follows (S,T) = with A being the 2-dimensional 



(3) 



describing the behavior of small perturbations b of the 
reference state. 

With respect to the magnetic energy balance Eq. (||) 
shows a remarkable difference to Eq. (||). Along with the 
term curl bx Bq which is again energy-conserving now as 
a second Hall-term curl Bo X b occurs which may well de- 
liver or consume energy (from/to b !) since in general the 
integral J v (curl -Bo x b) -curl b dV will not vanish. This re- 
flects the fact that the linearized Hall-induction equation 
describes the behavior of only a part of the total magnetic 
field. Actually, perturbations may grow only on expense 
of the energy stored in the background field. Considering 
(|J) , we can determine a scale below which the ohmic dis- 
sipation dominates the Hall-drift. Estimating | curl -Bo | 
and | curl b | by Bo and b/l, respectively, we find the crit- 
ical scale of b to be l cr u % 1/Bq, which is de-normalized 



lateral Laplacian (see p2[ ). With (|J) we obtain from (||) 
two coupled ordinary differential equations 



pt-t" +k 2 t 



= ik x f{s" - k 2 s) - ik x f"s 



ps — s" + k 2 s — ikyf's — —ik x ft , 



(6) 



L 



< 



L/(ujg T e ), identical with the expression derived 
in ||] considering the Hall-cascade in analogy with the 
turbulent flow of an incompressible fluid. 



where the dash denotes the derivative with respect to z. 
Together with appropriate boundary conditions Eqs. (||) 
define an eigenvalue problem with respect to p. 

We consider two types of boundary conditions (BC): 
For the vacuum condition we assume curl B = outside 
the slab and require continuity of all components of B 
across the boundary. For the perfect -conductor condi- 
tion an electric field must be prevented from penetrat- 
ing into the region outside the slab, that is, the normal 
magnetic and tangential electric field components must 
vanish at the boundary. In terms of the scalars s and t 
this means [s] = [s 1 ] = t = for the vacuum condition 
and s = t' = for the perfect conductor condition where 
[.] denotes the jump across a boundary. For t' = to 
be valid the vanishing of -Bq at the boundary is required. 
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Making use of the vacuum solutions vanishing at infin- 
ity for either halfspace, z > 1 and z < — 1, respectively, 
the vacuum boundary condition for s can be expressed 

as s' = ^fks for z = ±1, with k = \k\. 

Obviously, three distinguishable combinations of the 
boundary conditions are possible: vacuum on either side 
(VV), perfect conductor on either side (PP), vacuum on 
one and perfect conductor on the other side of the slab 
(PV). The latter choice comes closest to neutron star con- 
ditions if we think of the crust being neighboured upon a 
superconducting core on the one and a region with very 
low conductivity on the other side. The VV boundary 
condition may in turn be appropriate for a galactic disc. 
Since both the PV and the VV BC were to be consid- 
ered, we choose sufficiently curved background field pro- 
files, which obey them, i.e. f{z) = B (l +z)(l — z 2 ) and 
f(z) = B Q (l-z 2 ) for BC=PV and BC=VV , respectively. 

For certain ranges of the wave numbers k x , k y and for 
Bo > 3 we found eigenvalues p with a positive real part, 
i.e. exponentially growing perturbations. The depen- 
dence of the growth rate on the wave numbers for 
B = 1000 and BC = PV is shown in Fig. @. 

Figure |^ shows the dependence of growth rate and wave 
number k x of the fastest growing mode on Bq. 

An interesting feature is, that the maximum growth 
rates occur for all Bq considered at k y = 0. Of course 
this asymmetry is due to the choice of the background 
field: once it was chosen parallel to the y-direction the 
maximum growth rates would occur at k x — 0. Another 
interesting result is the dependence of the growth rate on 
the boundary conditions. BC = PV yields the largest val- 
ues, by a factor 1.6 ... 3 larger than for BC — VV, while 
BC = PP results in very small growth rates. Note that 
the most unstable eigenmodes are always non-oscillatory, 
though oscillating unstable ones exist. 

Evidently, the obtained growth rates are in agreement 
with the constraint, formulated above: In comparison 
with the background field decay the growth of the most 
unstable perturbations is a fast process; thus we may 
consider it as 'episodically unstable'. 

With respect to the asymptotic behavior a — > oo for a 
fixed (unnormalized!) background field one has to note 
that the time increment p is normalized on the ohmic 
decay rate (oc c -1 )- From Figure |^ it can be inferred 
9?(p) oc Bq , q < 1 for Bo > 100, which means that in the 
limit of negligible dissipation the growth rate in physical 
units tends to zero. 

Figure || shows the eigensolutions (s,t)(z) of the 
fastest growing mode for three different values of Bq and 
BC = PV. One can observe that with increasing Bq the 
toroidal field becomes more and more small-scaled and 
concentrated towards the vacuum boundary. In contrast, 
the corresponding poloidal field remains large-scaled. 

The magnetic field structure of the fastest growing 
mode for Bq = 2000 and BC = PV is shown in Fig. |. 

Clearly, any assignment of the results gained by help 
of a very simplified model to astrophysical objects has to 



be done with great care. Even when accepting the plane 
layer as a reasonable approximation of a neutron star's 
crust one has to concede that the very specific profiles of 
Bq assumed above may only exemplify the field structure 
in the crust. 

An acceptable approximation of the radial profile of a 
dipolar crustal field as given e.g. in jl3| will in general 
have to allow for a linear part and non-zero values at 
the boundaries. Moreover, the strong dependence of the 
conductive properties on the radial co-ordiante should 
anyway be taken into account. 

To get an impression of possible consequences for the 
evolution of neutron stars we now simply assume, that 
the real B-profile is sufficiently "curved" (i.e. its second 
derivative is big enough) and associate the parameter Bq 
with a typical value of the field. 

Assuming further electric conductivity and chemical 
composition to be constant, a = 5 x 10 26 s _1 and the 
relative atomic weight A/Z — 25, respectively, we find 
the normalization field at a density p = 10 14 g cm~ 3 to 
be 7 x 10 10 G (see e.g. pfl). That is, for typical (inner) 
crustal magnetic fields ranging between 7 x 10 12 G and 
1.4 x 10 14 G we find a Bq between 100 and 2000 and the e- 
folding time of the most rapidly growing unstable mode 
to be 0.0035 and 0.0003 times the Ohmic decay time, 
respectively. Thus, an initial perturbation will quickly 
evolve to a level at which the linear analysis is no longer 
feasible, that is, at which it starts to drain a remarkable 
amount of energy out of the background field. 

We want to emphasize again that a sufficient curvature 
of the background field profile is a necessary condition for 
the occurrence of an unstable behavior. Therefore nei- 
ther the derivation of the well-known helicoidal waves 
(whistlers) nor its modification presented in || could re- 
veal it because a homogeneous background field is used. 

With even more care we may speculate about possi- 
ble observational consequences. The instability discussed 
here may perhaps act effectively in the crust of not too 
young (age t ^ 10 5 yrs) neutron stars. For those stars, 
the small-scale field modes initially generated or existing 
in the crust have already been decayed and the magnetic 
field is concentrated almost completely in the large-scale, 
say, dipolar mode. Simultaneously, in the process of cool- 
ing the coefficient m*c/er e becomes smaller and smaller 
until the nonlinear Hall term in (|J) dominates the lin- 
ear ohmic term. From that moment the Hall-instability 
may rise small-scale modes down to scale lengths > l C rit 
on expense of the dipolar mode. This would lead to a 
change of the spin-down behaviour of isolated neutron 
stars. Deviations from the "standard" rotational evolu- 
tion will occur when the dipolar field decreases rapidly 
due to the instability. This may lead to the observation of 
braking indices n > 3 during the action of that insta- 
bility. Another possible observational consequence is due 
to enhanced Joule heating, which will keep the neutron 
star warmer than standard cooling calculations predict 
after an age critical for the onset of the Hall-instability 
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Qt 10 5 yrs). Third, the strong small-scale field compo- 
nents cause strong small-scale Lorentz forces which may 
be able to crack the crust. This could be observable in 
glitches or, depending on the available energy even in 
Gamma- and X-ray bursts |l5| . 

The authors are grateful to K.-H. Radler, D. Konenkov 
and K. Baumgartel (AIP), K. Gartner (Weierstrass Inst. 
Berlin), A. Reisenegger (PUC Chile) and D. Yakovlev 
(Joffe Inst. St. Petersburg) for helpful discussions. 



FIG. 1. Growth rate as a function of k x and k y for 
Bo = 1000 and BC = PV. Negative values were set to zero. 

FIG. 2. Growth rate and wave number fc™ ax of the fastest 
growing mode as functions of Bo. Solid and dashed lines 
correspond to BC = PV and BC = VV, while thick and thin 
lines correspond to growth rates and k x , respectively. 



FIG. 3. Moduli of (s,t)(z) of the fastest growing mode; 
BC = PV. Solid, dash-dotted and dashed lines refer to 
Bo = 2000, B = 100 and B = 10, respectively. 
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FIG. 4. Field structure of the fastest growing mode for 
Bo = 2000 and BC = PV. Arrows: b x ,z, grey shading: value 
of by, dark — into, light — out of the plane. 
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